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ABSTRACT 

We present a study of multiwavelength X-ray and weak lensing scaling relations for a sample of 
50 clusters of galaxies. Our analysis combines Chandra and XMM-Newton data using an energy- 
dependent cross-calibration. After considering a number of scaling relations, we find that gas mass 
is the most robust estimator of weak lensing mass, yielding 15 ± 6% intrinsic scatter at r^Q. The 
scatter does not change when measured within a fixed physical radius of 1 Mpc. Clusters with small 
BCG to X-ray peak offsets constitute a very regular population whose members have the same gas 
mass fractions and whose even smaller (< 10%) deviations from regularity can be ascribed to line 
of sight geometrical effects alone. Cool-core clusters, while a somewhat different population, also 
show the same (< 10%) scatter in the gas mass-lensing mass relation. There is a good correlation 
and a hint of bimodality in the plane defined by BCG offset and central entropy (or central cooling 
time). The pseudo-pressure Yx does not discriminate between the more relaxed and less relaxed 
populations, making it perhaps the more even-handed mass proxy for surveys. Overall, hydrostatic 
masses underestimate weak lensing masses by 10% on the average at r^^; but cool-core clusters are 
consistent with no bias, while non-cool-core clusters have a large and constant 15 — 20% bias between 
r^QQ and r^Q, in agreement with N-body simulations incorporating unthermalized gas. For non-cool- 
core clusters, the bias correlates well with BCG ellipticity. We also examine centroid shift variance 
and and power ratios to quantify substructure; these quantities do not correlate with residuals in the 
scaling relations. Individual clusters have for the most part forgotten the source of their departures 
from self-similarity. 

Subject headings: Dark matter — X-rays: galaxies: clusters — gravitational lensing 



1. INTRODUCTION 

Within the context of the currently favored hierarchi- 
cal model for structure formation, massive clusters of 
galaxies are, as a population, the most recently formed 
gravitationally bound structures in the cosmos. Conse- 
quently, characteristics such as the shape and evolution- 
ary behavior of their mass function can, in principle, be 
exploited as precision probes of cosmology. The resulting 
estimates of parameters — such as the amplitude of the 
primordial fluctuations and the density and equation of 
state of the mysterious dark energy — can certainly com- 
plement and even compete with determinations based on 
stu dies of t he cosmic microwave background (for a review 
see 



Allen et al. 201 lai) 



The etticacy oi clusters as cosmological probes depends 
on three factors: (1) the ability to compile a large well- 
understood catalog of clusters; (2) the identification of 



an easily determined survey observable (or combinations 
thereof) — hereafter referred to as a "mass proxy" — 
that can offer an accurate measure of cluster masses; and 
(3) the existence of a well-calibrated relationship between 
the mass proxy and the actual mass of the cluster. Of 
these, we shall focus our attention on the latter two since 
at present, the effective use of clusters as cosmological 
probes is primarily limited by systematic errors in the 
esti mates of the true mass of the cluster (He nry et al. 
"2009 ; Vikh linin et al.||2009bj |Mantz et al.||2010[ ). 

One of the first — and still among the rnost commonly 
used — mass proxies is the "hydrostatic mass estimate", 
derived from X-ray observations under the assumption 
that the clusters are spherically symmetric and that the 
hot, diffuse. X-ray emitting gas in galaxy clusters is 
in thermal pressure-supported hydrostatic equilibrium 
(HSE). Over the years, mismatch between hydrostatic 
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mass estimates and mass estimates derived by alternate 
means have led a num ber of researchers to question the 
use of this proxy (e .g Miralda-Escude & Babul"1995', 'Fis- 
cher fc Tyson|1997i Girardi ct al. 1997; Ota ct al. 2004). 
Kecent studies suggest that the HSE masses of relaxed 
clusters are subject to a systematic 10%-20% underesti- 
mate which grows to 30% or more for unrelaxed system s 



( Arnaud et al.|2007[ [Mahdavi et al.|2008[|Lau et al.|2009[ ). 
iNumerical siinulation studies suggest that this bias is due 
to incomplete thcrmalization of the hot diffuse intraclus- 
ter medium (ICM) (Evrard 1990; Nagai et al.|2007[|Shaw 
et al.||2010| (Rasia et al. 2Ul'2f~' 

Concerns with the HSE mass estimate have renewed in- 
terest in identifying more well-behaved mass proxies that 
can give unbiased estimates of the cluster mass. One ex- 
ample of such an X-ray mass proxy is Yx, the product 
of the gas mas s M„ and ICM te mperature Tx within a 
given aperture (Kravtsov et al. 2006) . In numerical simu- 
lation studies, this pressure- like quantity has been shown 
be a much better mass proxy and has been successfully 
deployed in measurements of cosmological p arameters 
including the dark energy equation of state ( Vikhlinin 
et al.||2009a|b I . More recently, the gas mass Mg has also 
emerge d as a mass proxy with similar predict ive power 



to Yx (Okabe et al. 2010 



Allen et al. 2011b). Success 



in tests involving simulated clusters is necessary but far 
from sufficient. At present, numerically simulated clus- 
ters capture only a fraction of the physical processes that 
affect the intracluster medium in real clusters. 

An alternative way of independently testing the valid- 
ity of the individual mass proxies is via multiwavelength 
observations. Specifically, comparisons of X-ray proxies 
and weak gravitational lensing masses (M^) are partic- 
ularly interesting given the fact that gravitational lens- 
ing provides a total mass estimate that neither depends 
on baryonic physics nor requires any strong assumptions 
about the equilibrium state of the gas and dark matter, 
and which can be determined over a wide range of spa- 
tial scales. However, lensing measures the projected (2D) 
mass and converting this to a unprojected (3D) mass is 
subject to systematic and statistical errors related to the 
geometry of the mass distribution, its orientation along 
the line of sight, and projection of extra-cluster mass 
along the line of sight. In individual cases, these effects 
can result in an under- or over- estimate of the cluste r 
mass of as much as a factor of 2 ( Feroz &: Hobson||2012 ), 
depending on the specific technique used. However, nu- 
merical simulations find that averaging over a sample of 
clusters, 3D weak lensing m ass estimates are typically 
shghtly biased low (^3-6%) ([Becker k Kravt sov 2011] 



Bahe et al.][2C)12[ [High et al [|2012[ [KasiT^ al.j|2012| 
where the actual amount depends on the range of phys- 
ical radii used in the weak lensing analysis. 

As such, weak lensing masses are, at present, the best 
measures of cluster mass and very well suited for use in 
calibrating the different mass proxies and identifying the 
best one of the lot. Moreover, the study of the rela- 
tionship between the weak lensing mass estimate and an 
observable mass proxy can potentially yield important 
insights into the physics at play within cluster environ- 
ments. These are the goals of the present paper. 

To facilitate our study, we have assembled a sample of 
galaxy clusters named the Canadian Cluster Comparison 



Projeclj^ We describe this sample in sj2] In the present 
study, we restrict ourselves to studying the relationships 
between weak lensing mass determinations and the mass 
proxies derived jointly from Chandra and XMM-Newton 
observations. We use the Jo int Analysis o f Observations 
(JACO) code base (Mahdavi et al.||2007[ ) to derive the 
mass proxies of interest from the X-ray data. JACO 
makes maximal use of the available data while incorpo- 
rating detailed corrections for instrumental effects (for 
example, we model spatial and energy variations of the 
PSF for both Chandra and XMM-Newton) to yield self- 
consistent radial profiles for both the dark and the bary- 
onic components. Further details are given in §2.4[ In ^ 



we summarize our data reduction procedure; in ^ 2.4 



we 



describe our mass modeling technique. Our quantitative 
measures of substructure, the luminosity-temperature re- 
lation, the lensing mass-observable relations, and devia- 
tions from hydrostatic equilibrium are discussed in < 



j]5] and Sj6j respectively. We conclude in Throng' 
bnc paper we take — 70 km/s/Mpc, TFm = 0.3, and 



Oa = 0.7. 

2. SAMPLE AND DATA REDUCTION 
2.1. Sample Characterization 

The Canadian Cluster Comparison Project (CCCP) 
was established primarily to study the different baryonic 
tracers of cluster mass and to explore insights about the 
thermal properties of the hot diffuse gas and the dynami- 
cal states of the clusters that can be gained from cluster- 
to-cluster variations in these relationships. 

For this purpose, we assembled a sample of 50 clusters 
of galaxies in the redshift range 0.15 < z < 0.55. Since 
we wanted to carry out a weak lensing analysis, we re- 
quired that the clusters be observable from the Canada- 
France-Hawaii Telescope (CFHT) so we could take ad- 
vantage of the excellent capabilities of this facility. The 
latter constraint restricts our cluster sample to systems 
at —15° < declination < 65°. We also required our 
clusters to have an ASCA temperature ksTx > 3 keV. 
To establish cluster temperature, we pri marily relied on 
a systematically reduced cluster catalog of Horner ( 2001 1 
based on ASCA archival data, although in a tew instances 
we used temperatures from other (published) sources. 

As a starting point, we scoured the CFHT archives for 
clusters with high quality optical data suitable for weak 
lensing analysis, including observations in two bands. We 
identified 20 suitable clusters observed with the CFH12k 
camera and with B and R band data meeting our cri- 
teria. Nearly half of these clusters were originally ob- 
served as part of the Canadian Network for Observational 
Cosm ology (CNOCl) Survey ([Yee et al.[[1996| [Carlberg 



et al. 1996) and comprise the brightest clusters in the 
em Obs ervatory Extended Medium Sensitivity Sur 



hlinst 

vey (EMSS) ( [Gioia et al. 

pie is known to have a mi ^ , 

clusters wit h pronounced substructure ( Pesce et al.|1990 



1990[). Since the EMSS sam- 
bias agains t X-ray luminous 



Donahue et al. 1992 ; Ebcling et al. 2000p , and we were 



specifically interested m putting together a representa- 
tive sample of clusters that encompassed the spectrum of 
observed variations in thermal and dynamical states, we 
randomly selected 30 additional clusters from the Horner 

^ Not to be confused with the Chandra Cluster Cosmology 
Project, which forms an identical acronym 
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sample that met our temperature, declination and red- 
shift constraints and additionally, guaranteed that our 
final sample fully sampled the scatter in the Lx vs. Tx 
plane. Of these systems, those without deep, high qual- 
ity optical data were observed with the CFHT MegaCam 
wide-field imager, using the g' and r' optical filter sets. 
The resulti ng weak lensing masse s for this sample are 



discussed in Hoekstra et al. (2012 1. 

Our final sample comprises 5U clusters listed in Table 1. 
All except 3 clusters have been observed by the Chandra 
Observatory. These three, plus 21 others, have also been 



sample have been used in several prior studies ( Hoekstra 


2007 


Mahdavi et al. 


2008 


Bildfell et al. 


200Br 


2012). 



i'he CCCP sample has served as the source for studies of 
individual clusters that are interesting in t heir own right, 
such asAbcU 520 and IRAS 09104-^4 109 (Mahdavi et al. 
2007||Jee et al. 2012; O'SuUivan et al..2012j 



Tn the left panel of t'igure 1, we conipare the distribu- 
tion of the CCCP clusters in the Lx — Tx plane to those 
of two better characterized samples of galaxies clusters: 
MACS ([Eb cling et al.|[20To|) and HIFLUGCS ( [Reiprich 
fc Bohrmg cr 2002), both of which employ well-dehned 
Ilux- based selection criteria based on the ROSAT All- 
Sky Survey. HIFLUGS is on the average a lower redshift 
sample compared to our CCCP sample, and MACS is 
on the average at a higher redshift. The samples have 
comparable scatter, suggesting that our CCCP sample is 
not significantly more biased than HIFLUGCS or MACS, 
which have better understood selection functions. In 
the right panel of Figure 1, we plot the distribution of 
the orthogonal scatter about the mean Lx-Tx of the 
all three samples combined. A KS test indicates that 
the three distributions are statistically indistinguishable. 
This confirms that while the CCCP sample may not be 
a complete sample, it is a representative sample in that 
it properly captures the scatter in the Lx — Tx and to 
the extent that these have physical origins, the range of 
cluster thermal and dynamical states. 

2.2. Choice of density contrast 

For most of what follows, we study masses, tempera- 
tures, substructure measures, and other thermodynamic 
quantities integrated within a specific spherical radius. 
The choice of this radius is not obvious; using fixed physi- 
cal radii has the advantage of straighforwardness, but the 
disadvantage that we would be probing characteristically 
different regions of clusters as a function of masses. Using 
fixed overdensity radii ta (defined such that ta contains 
a mean matter density of A times the critical density 
of the universe at the redshift of the cluster) is a bet- 
ter choice, but even here, the value of A to use is not 
quite obvious. At the redshift of our sample, X-ray data 
quality tends to be best around ^2500, but most of the 
literature lists properties at r5oo- Even after a choice of 
A, one must still decide whether to use the lensing or 
X-ray value, since they are not guaranteed to agree. 

We choose to standardize the bulk of our discussion on 
the weak- lensing overdensity radius r^Q, because lensing 
masses a re likely to be more un biased for non-relaxed 



clusters (Meneghetti et al. 2010). For the most part 



our results do not signihcantly change if we switch to X 
ray ^500; one exception is the mass-temperature relation 
below, which tightens significantly with the switch. In 



|5.3[ we also consider scaling relations with observables 
measured within fixed physical radii, because these are 
more likely to be useful for calibrating large data sets. 

2.3. Weak Lensing Overview 



The cluste rs in our sample were drawn from |Hoekstra| 
(2012), which contains a weak lensing analysis of 



2k and Megacam data from the Canada-France- 
Hawaii Telescope . We refer interested readers to HoeF] 
stra et al. (2012) for details of the data reduction and 
weak lensing analysis procedure. 

We base our lensing masses on the aperture mass esti 
mates (for details see the discussion in 
2007]) 
b 



3.5 in Hoekstra 

This approach has the advantage that it is prac- 
tically model independent. Additionally, as the mass es- 
timate relies only on shear measurements at l arge radii 



contam inat ion by cluster member s is minimal. Hoekstra 
(2007) and Hoekstra et al. (2012) removed galaxies that 



he on the cluster red-sequence and boosted the signal 
based on excess number counts of galaxies. As an ex- 
treme scenario we omitted those corrections and found 
that the len sing masses change by only a few percent; for 
details see ( Hoekstra et al.|[2012 ). Hence our masses are 
robust against contamination by cluster members at the 
percent level. 

The weak lensing signal, however, only provides a di- 
rect estimate of the projected mass. To calculate 3D 
masses from the model-independent 2D aperture masses 
we project and renormalize a densitv profile of th e form 



Ptot{r) oc r"i(r2oo + cr)"^ (Navarro et al. 1997). The 



relationship between the concentration c and the virial 
mass is fixed at c oc MZ^^^^{1 + z) from numerical sim- 



ulations (Duffy et al. 2008). Hence, the deprojection 
itself, though well motivated based on numerical sim- 
ulations, is model dependent. However, the model de- 
pendence is weak — 20% variations in the normalization 
of the mass-concentration relationshi p yield » 5% varia- 
tions in the measured masses ( .Hoekstra et al.|2012[ §4.3). 
We also note that the lensing analysis ditters from the 
X-ray analysis in that in the X-ray analysis, no mass- 
concentration relationship is assumed (i.e., the concen- 
trations and masses are allowed to vary independently). 
We plan to address the effects of relaxing the lensing 
mass-concentration relation in a future paper. 

2.4. X-ray Data Reduction 



We refer the reader to Mahdavi et al. (2007) for de- 
tails of the X-ray data reduction procedure, which we 
briefly summarize and update here. We use both Chan- 
dra CALDB 4.2.2 (April 2010) and CALDB 4.4.7 (De- 
cember 2011). We also check our results against the lat- 
est CALDB (4.5.1) at the time of writing. For XMM- 
Newton we use calibration files up-to-date to January 
2012; we also checked calibration files dating as far back 
as April 2010. We detected no statistically significant 
changes in the calibration files over this period fo r cit her 
Chandra or XMM-Newton, except as detailed in \ 2.6 be- 
low. 

We follow a standard data reduction procedure. We 
use the software packages CIAO (Chandra) and SAS 
(XMM-Newton) to process raw event files using the rec- 
ommended settings for each observation mode and de- 
tector temperature. Where possible, we make event 
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Tjj (keV) Orthogonal Dist. (dex) 



Fig. 1. — Comparison of the luminosity-temperature relationship for JACO/CCCP sample (solid dots), HIFLUGS (open dots) and MACS 
(stars) 




Fig. 2. — (top) (ie/t)Comparison of XMM-Newton and Chandra X-ray masses within lensing J'^qq. (right) When scaling the Chandra 
effective area by a power law in energy of slope <J = 0.07, Chandra and XMM-Newton masses come into agreement. 



JACO/CCCP Clusters of Galaxies 
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grade selections that maximize the data quality for ex- 
tended sources (including the VFAINT mode optimiza- 
tions for Chandra). We use the wavelet detection al- 
gorithm WAVDETECT on exposure-corrected images to 
identify contaminating sources; we masked out point and 
extended sources using the detected wavelet radius. Each 
masking was checked by eye for missing extended sources 
or underestimated masking radii. 

The bulk of the X-ray background consists of a parti- 
cle component which bypasses the mirror assembly, plus 
an astrophysical component that is folded through the 
mirror response. To remove the particle background we 
match the 8-12 keV photon count rate from the outer re- 
gions of each detector to the recommended blank sky ob- 
servations for each detector, and then subtract the renor- 
malized blank-sky spectra. What remains is the source 
plus an over- or under-subtracted astrophysical back- 
ground, plus in some cases residual particle background. 
All these residual backgrounds are modeled jointly with 
the spatially resolved ICM model spectra, and their pa- 
rameters marginalized over for the final results. 

To extract spatially resolved spectra, we find the sur- 
face brightness peak in the Chandra image (if available) 
or XMM-Newton image (if Chandra is not available). 
We then draw circular annuli that contain a minimum of 
1500 background-subtracted photon counts; where both 
Chandra and XMM-Newton data are available, the an- 
nuli are taken to be exactly the same for both sets of 
observations, with the minimum count requirement be- 
ing imposed on the Chandra data (for photons within 
8') or XMM-Newton data (for photons outside 8')- We 
then compute appropriately weighted ancilliary response 
files (ARE) and redistribution matrix files (RME) for 
each spectrum, and subtract appropriately scaled parti- 
cle background spectra. We emphasize that all spectra 
for each cluster undergo a simultaneous joint fit using a 
forward-convolved spectral model of the entire cluster, so 
that the choice of 1500 background-subtracted counts is 
not a sensitivity-limiting factor. That is to say, in no case 
is a single measurement derived from a single spectrum 
of 1500 counts, but rather such spectra are fit together 
in large batches on a cluster-by-cluster basis. 

The detailed properties of the sample, including global 
X-ray temperatures and bolometric X-ray luminosities, 
masses, and substructure measures are listed in tables [T] 
and El 

2.5. X-ray Mass Modeling 

Here we summarize and upd ate the modeling proce- 
dure of Mahdavi et al. (2007), in which the cluster is 
spherically symmetric and that the gas is in thermal 
pressure supported hydrostatic equilibrium within the 
cluster potential. The essence of the technique is to di- 
rectly compare the observed spatially resolved spectra 
with model predictions. Eor a spectrum observed in an 
annulus with inner and outer radii Ri and i?2 , the model 



R2 



2TrRdR 



UennK [T[r), Z{r)] 



2rdr 



- i?2 

where r denotes unprojected radius, R denotes projected 
radius, Tmax is the termination radius of the X-ray gas 
(taken to be rioo in this paper), A^ is the frequency- 



dependent cooling function which is a function of tem- 
perature T and metallicity Z, and rie and uh are the 
electron and hydrogen number density, respectively. 

One feature of the above method is that the unpro- 
jected temperature profile is calculated self-consistently 
assuming hydrostatic equilibrium of assumed gas and 
dark matter density profiles. As a result, we never have 
to specify or fit a temperature profile; temperature is 
merely an intermediate "dummy" quantity connecting 
the gas and dark matter mass distributions to the X-ray 
spectra. This avoids subjective weighting involve d in the 
fitting of 2p projected temperature profiles ( |Mazzotta 



et al. 2004 Rasia et al. 2005 



Vikhlinin|2006 ), which are 



more ditticult to correct tor the ettects of PSl' distortion 

2.6. Parameters of the Hydrostatic Model 

The hydrostatic model assumes a flexible spherical 
electron density distribution 



n{r) : 



S(r,r,„,/3o) + 



ne,B{r, r^^,l3i) + ne^B{r, r^^,l32) 
where the familiar "beta" model is 



(2) 
(3) 

(4) 



In other words, the gas mass profile consists of a fully 
general triple "beta" model profile, where the first beta 
model is further allowed to be multiplied by a single 
pow er law. The metallici ty distribution is modeled as 
(e.g. [Pizzolato et"al]|2003[ ) 

with r^, Pzi and Zq free parameters. Finally, the total 
mas s distribution (baryon s and dark matter) are modeled 
as a 



Navarro et al.| (|1997|) profile: 

Mo 



Ptot 



r{cr -|- ta)^ 



(6) 



where Mq is the normalization, c is the halo concentra- 
tion, and TA is the overdensity radius (see above). These 
are also free parameters, except that rather than fitting 
Mq, we fit Ma — the mass within ta — as the normaliza- 
tion constant (because there is a one-to-one relationship 
between Mq and Ma). 

2.7. Joint Calibration of Chandra and XMM-Newton 

Masses 

Where available, we use both Chandra and XMM- 
Newton data for a cluster. This has several advantages: 
in the inner regions, Chandra is able to resolve the clus- 
ter cores well; while XMM-Newton's wider field of view 
yields better coverage of the outer regions of the clus- 
ter. The simultaneous coverage of intermediate regions 
helps constrain residual backgrounds following blank sky 
subtraction. 

When combining Chandra and XMM-Newton data, 
cross-calibration is a significant issue. In general, there 
are slight differences among the responses of the Chan- 
dra ACIS and the XMM-Newton pn, MOSl, and M0S2 
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TABLE 1 
Basic Properties of the Sample 



Cluster 


RA 


DEC 


z 


Chandra 


Exposure 


XMM-Newton 


Exposure 


-^X.all bol 500 




Name 


J2000 


J2000 




ObsID 


s 


ObsID 


s 


10* erg s~^ 


ke'v 


3C295 


14:11:20.52 


-f52:12:09.9 


0.464 


2254 


87914 






1.77 ± 0.06 


5.9 ±0.6 


Abell0068 


00:37:06.65 


-1-09:09:24.0 


0.255 


3250 


9986 


0084230201 


14068 


1.87 ± 0.05 


6.8 ±0.4 


AbellOllSN 


00:55:50.37 


-f 26:24:36. 6 


0.197 


3233 


49719 


0203220101 


21393 


1.00 ±0.01 


5.2 ±0.1 


AbellOllSS 


00:56:00.17 


-f 26:20:29. 5 


0.197 


3233 


49719 


0203220101 


21309 


0.68 ±0.02 


5.5 ±0.3 


Abell0209 


01:31:53.42 


-13:36:46.3 


0.206 


3579 


9986 


0084230301 


11219 


1.86 ±0.03 


7.0 ±0.3 


Abell0222 


01:37:34.25 


-12:59:30.8 


0.207 


4967 


45078 


0502020201 


23178 


0.50 ±0.02 


4.1 ± 0.3 


Abell0223S 


01:37:56.06 


-12:49:12.8 


0.207 


4967 


45078 


0502020201 


23206 


0.48 ±0.01 


5.6 ±0.3 


Abell0267 


01:52:42.38 


-1-01:00:48.0 


0.231 


3580 


19624 


0084230401 


10421 


1.47 ± 0.04 


6.8 ±0.3 


Abell0370 


02:39:53.18 


-01:34:34.9 


0.375 


515 


68532 






1.89 ±0.05 


7.4 ± 0.6 


Abell0383 


02:48:03.33 


-03:31:45.1 


0.187 


2320 


19285 


0084230501 


20237 


1.51 ± 0.01 


3.9 ±0.1 


Abell0520 


04:54:10.10 


-1-02:55:18.3 


0.199 


4215 


66274 


0201510101 


21915 


1.75 ±0.04 


7.8 ±0.4 


Abell0521 


04:54:06.30 


-10:13:16.9 


0.253 


901 


38626 






1.09 ±0.03 


5.9 ±0.3 


Abell0586 


07:32:20.16 


-1-31:37:56.6 


0.171 


530 


10043 






1.62 ±0.06 


5.4 ± 0.4 


AbellOeil 


08:00:56.96 


-1-36:03:22.0 


0.288 


3194 


36114 






1.94 ± 0.06 


7.0 ±0.9 


Abell0697 


08:42:57.29 


-1-36:21:56.2 


0.282 


4217 


19516 






3.32 ±0.10 


10.0 ± 1.1 


Abell0851 


09:43:00.39 


-1-46:59:20.4 


0.407 






0106460101 


15731 


0.91 ± 0.03 


5.7 ± 0.5 


Abell0959 


10:17:35.61 


-f59:33:53.4 


0.286 






0406630201 


4134 


0.74 ± 0.04 


6.5 ± 1.7 


Abell0963 


10:17:03.63 


-1-39:02:48.3 


0.206 


903 


36289 


0084230701 


17234 


1.96 ±0.04 


6.2 ±0.2 


Abelll689 


13:11:29.52 


-01:20:29.8 


0.183 


6930 


76144 


0093030101 


24457 


4.91 ± 0.02 


9.1 ± 0.2 


Abelll758E 


13:32:46.43 


-1-50:32:25.9 


0.279 


2213 


55220 






1.70 ±0.03 


9.6 ±0.9 


Abelll758W 


13:32:38.70 


-f50:33:23.0 


0.279 


2213 


55220 






1.30 ±0.03 


9.9 ± 1.4 


Abelll763 


13:35:18.16 


-1-40:59:57.7 


0.223 


3591 


19595 


0084230901 


8852 


2.01 ± 0.04 


7.0 ±0.3 


Abelll835 


14:01:01.90 


-1-02:52:42.7 


0.253 


6880 


117918 


0098010101 


16021 


7.00 ±0.03 


7.0 ±0.1 


Abelll914 


14:26:02.80 


-1-37:49:27.3 


0.171 


3593 


18865 


0112230201 


17025 


3.61 ± 0.06 


9.3 ±0.3 


Abelll942 


14:38:21.90 


-1-03:40:12.9 


0.224 


3290 


55716 






0.44 ± 0.02 


4.4 ± 0.5 


Abell2104 


15:40:08.09 


-03:18:16.5 


0.153 


895 


49199 






1.62 ±0.02 


5.9 ±0.3 


Abell2111 


15:39:41.74 


-1-34:25:01.9 


0.229 


544 


10299 






1.13 ±0.04 


6.0 ± 1.0 


Abell2163 


16:15:46.05 


-06:09:02.6 


0.203 


1653 


71148 






5.95 ±0.10 


11.0 ±0.4 


Abell2204 


16:32:46.92 


-f05:34:32.4 


0.152 


7940 


77141 


0306490201 


13093 


5.34 ± 0.02 


7.1 ± 0.2 


Abell2218 


16:35:50.89 


-1-66:12:36.9 


0.176 


1666 


30693 


0112980101 


13111 


1.70 ±0.02 


6.8 ±0.2 


Abell2219 


16:40:20.20 


-1-46:42:35.3 


0.226 


896 


42295 






4.66 ±0.08 


8.9 ±0.6 


Abell2259 


17:20:07.75 


-f 27:40:14.7 


0.164 


3245 


9986 






1.00 ±0.07 


5.3 ±0.6 


Abell2261 


17:22:27.12 


-1-32:07:58.9 


0.224 


5007 


24316 






3.29 ±0.09 


6.4 ± 0.5 


Abell2390 


21:53:36.82 


-1-17:41:44.7 


0.228 


4193 


93782 


0111270101 


8100 


5.39 ±0.03 


8.8 ±0.2 


Abell2537 


23:08:22.23 


-02:11:30.3 


0.295 


4962 


36193 


0205330501 


6267 


1.78 ±0.06 


6.8 ±0.7 


CL0024.0+1652 


00:26:35.94 


-1-17:09:46.2 


0.390 


929 


39417 






0.49 ± 0.03 


4.6 ± 1.1 


MACSJ0717.5+3745 


07:17:31.39 


-1-37:45:24.8 


0.548 


4200 


58912 






6.08 ±0.15 


11. 3± 1.0 


MACSJ0913. 7+4056 


09:13:45.49 


-f 40:56:28. 7 


0.442 


10445 


76159 






3.18 ±0.04 


6.0 ±0.3 


MS0015.9+1609 


00:18:33.74 


-1-16:26:09.0 


0.541 


520 


67410 


0111000101 


22477 


4.22 ±0.08 


8.7 ± 0.7 


MS0440.5+0204 


04:43:09.99 


-1-02:10:19.3 


0.190 


4196 


22262 






0.59 ±0.03 


3.4 ± 0.4 


MS0451. 6-0305 


04:54:11.24 


-03:00:57.3 


0.550 


902 


43420 






3.96 ±0.12 


10.2 ± 1.2 


MS0906.5+1110 


09:09:12.73 


-1-10:58:28.4 


0.174 


924 


29752 






1.07 ± 0.04 


5.5 ±0.3 


MS1008. 1-1224 


10:10:32.52 


-12:39:53.1 


0.301 


926 


25222 






1.11 ± 0.04 


5.8 ±0.6 


MS1231.3-I-1542 


12:33:55.01 


-1-15:26:02.3 


0.233 






0404120101 


26520 


0.34 ± 0.01 


4.6 ±0.4 


MS1358. 1-1-6245 


13:59:50.56 


-f62:31:05.3 


0.328 


516 


50989 






1.69 ±0.04 


6.6 ±0.6 


MS1455.0-I-2232 


14:57:15.05 


-f 22:20:33. 2 


0.258 


4192 


91626 


0108670201 


22571 


3.14 ± 0.02 


4.4 ± 0.1 


MS1512.4-I-3647 


15:14:22.47 


-1-36:36:20.9 


0.372 


800 


36400 






0.78 ±0.03 


3.1 ± 0.3 


MS1621.5-I-2640 


16:23:35.05 


-1-26:34:22.1 


0.426 


546 


30062 






1.03 ±0.07 


6.5 ± 1.4 


RXJ1347.5-1145 


13:47:30.59 


-11:45:09.8 


0.451 


3592 


57458 


0112960101 


21712 


13.15 ±0.22 


12.1 ± 0.4 


RXJ1524. 6-1-0957 


15:24:38.85 


-1-09:57:41.8 


0.520 


1664 


49849 






0.48 ± 0.04 


5.1 ± 1.0 



Note. — i'X.all,bol.500 the bolometric X-ray luminosity and Taii.500 is the temperature measured using aU data within )"Jq()" of the 
cluster center. 



detectors. Even after over a decade in flight, the source 
of these differences has not been conclusively identified. 
Typically, compariso ns show that Chandra t emperatures 



are 5 — 15% higher (Snowden et al. 2008 Reese et al. 



2010|. The most recent calibration tests ( [I'sujimoto 



2011 1 use the G21.5-0.9 pulsar (which is fainter 
the Crab nebula, and hence not 



et al 

than the usual source 

subject to detector pileup). [Tsujimoto et al.| (2011| find 
that the XMM-Newton pn has a 157o lower nux in the 
2.0-8.0 keV energy band compared to t he Chandra ACIS- 
S. Th is confirms an earlier finding by [Nevalainen et ah] 
(2010) who found similar results. Lower hard band tiux 



naturally leads to lower X-ray temperatures when 0.5- 
2.0 keV photons are also included. This primarily affects 
masses for which spectral line emission is not dominant 



(i.e., in hot, k T > 4 keV clusters). It is at this point 
unknown where the source of the disagreement lies and 
which instrument is better calibrated. 

Figure [2] shows the X-ray mass measured within lens- 
ing r^QQ for the 19 clusters in our sample which contain 
both Chandra and XMM-Newton data. Shown are the 
resuhs for CALDB 4.2.2 (April 2010). We also checked 
CALDB 4.4.7 (December 2011) and CALDB 4.5.1 (June 
2012). The calibration for our sample changed little 
during this period, and in all three cases, we find that 
Chandra masses are higher than XMM-Newton masses 
by roughly 15%. All observations were recorded prior 
to 2010, and taken as a whole, the change in the Chan- 
dra masses of these systems is not statistically significant 
between CALDB 4.2.2 and 4.4.7. We adopt the 2010 
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TABLE 2 

Mass and Substructure Properties at rsoo 



Cluster 




' 500 








'^^hydro 
1O"M0 




Dbcg 


wx 


P3/P0 






Mpc 




10"Mq 


keV cm^ 


kpc 


103 X 


xlO"'^ 


3C295 


1 


.06 ± 0.06 


5.7+ 1, 


.2 


0.62 + 0.03 


3.9 + 1.0 


2.6 + 0.8 


12 


5.7 + 1.4 


0.25 + 0.20 


Abell0068 


1 


.16 ± 0.09 


5.9+ 1, 


.6 


0.77 + 0.01 


5.1 + 1.0 


35.1 + 5.8 


15 


14.0 + 2.1 


0.74 + 0.60 


Abell0115N 


1 


.03 ± 0.12 


3.9+ 1, 


.5 


0.60 + 0.01 


4.1 + 0.2 


4.7 + 0.7 


10 


59.6 + 0.8 


3.31 + 0.82 


AbellOllSS 


1 


.14 + 0.07 


5.3+ 1, 


.2 


0.80 + 0.01 


4.2 + 0.3 


36.9 + 9.8 


143 






Abell0209 


1 


.24 + 0.07 


6.8+ 1, 


.4 


1.02 + 0.02 


5.6 + 1.1 


21.7 + 3.2 


16 


7.1 + 1.1 


1.47 + 1.10 


Abell0222 


1 


.16 + 0.07 


5.7+ 1, 


.3 


0.61 + 0.01 


2.4 + 0.6 


56.9 + 9.7 


10 


38.4 + 3.2 


1.21 + 0.96 


Abell0223S 


1 


.24 + 0.10 


6.9 + 2, 


.0 


0.68 + 0.04 


3.3 + 1.6 


19.0 + 3.3 


8 


36.9 + 2.4 


2.26 + 1.62 


Abell0267 


1 


.13 + 0.09 


5.3+ 1, 


.5 


0.63 + 0.01 


5.7 + 0.6 


25.2 + 3.8 


77 


22.0 + 1.7 


0.27 + 0.22 


Abell0370 


1 


.43 ± 0.06 


12.8+2, 


.0 


1.00 + 0.05 


8.6 + 6.0 


77.1 + 26.2 


23 


15.8 + 1.3 


0.62 + 0.49 


AbellOSSS 


1 


.04 + 0.13 


4.0+ 1, 


.8 


0.39 + 0.01 


4.6 + 0.6 


6.3 + 0.8 


< 3 


3.1 + 0.6 


0.36 + 0.26 


Abell0520 


1 


.16 + 0.07 


5.6+ 1, 


.3 


0.85 + 0.01 


7.3 + 0.3 


82.8 + 6.5 


341 


100.7 + 0.7 


4.66 + 0.97 


Abell0521 


1 


.19 + 0.08 


6.3+ 1, 


.6 


1.06 + 0.02 


5.0 + 1.3 


12.8 + 3.4 


33 


58.5 + 1.7 


8.59 + 3.41 


Abell0586 


1 


.18 + 0.09 


5.6+ 1, 


.6 


0.65 + 0.08 


3.9 + 0.6 


31.1 + 8.7 


11 


7.3 + 1.6 


0.59 + 0.49 


AbcllOeil 


1 


.13 + 0.06 


5.7+ 1, 


.3 


0.66 + 0.05 


6.0 + 0.9 


5.0 + 2.0 


4 


8.0 + 0.8 


0.68 + 0.42 


Abell0697 


1 


.35 + 0.05 


9.7+ 1, 


.3 


1.56 + 0.03 


10.9 + 1.5 


22.6 + 6.5 


20 


9.1 + 1.1 


0.20 + 0.16 


Abell0851 


1 


.32 + 0.09 


10.5 + 2, 


.5 


0.97 + 0.02 


7.4 + 2.3 


73.0+ 18.2 


278 


30.5 + 3.4 


13.15 + 7.51 


Abell0959 


1 


.26 + 0.07 


7.8+ 1, 


.7 


0.75 + 0.03 


5.6 + 0.5 


34.9 + 11.4 


36 


42.6 + 4.5 


7.70 + 6.55 


Abell0963 


1 


.00 + 0.10 


3.7+ 1, 


.3 


0.57 + 0.01 


4.7 + 0.5 


10.4 + 1.1 


6 


5.5 + 0.5 


0.12 + 0.10 


Abelll689 


1 


.57 + 0.09 


13.7+2, 


.7 


1.27 + 0.01 


9.7 + 0.6 


8.2 + 1.0 


5 


4.1 + 0.3 


0.08 + 0.04 


Abelll758E 


1 


.37 + 0.08 


10.1 + 2, 


.3 


1.23 + 0.04 


9.4 + 0.6 


20.1 + 6.9 


25 


117.9 + 1.2 


8.42 + 1.62 


Abelll758W 


1 


.37 + 0.06 


10.0+ 1, 


.4 


0.93 + 0.07 


11.5 + 1.6 


21.1 + 6.2 


25 


117.9 + 1.2 


8.42 + 1.62 


Abelll763 


1 


.40 + 0.10 


10.1 + 2, 


.5 


1.34 + 0.01 


3.9 + 0.7 


70.0+ 13.8 


7 


22.9 + 1.2 


0.97 + 0.64 


Abelll835 


1 


.30 + 0.05 


8.4+ 1, 


.3 


1.21 + 0.01 


9.9 + 0.7 


2.2 + 0.1 


6 


3.9 + 0.2 


< 0.1 


Abelll914 


1 


.18 + 0.05 


5.6+ 1, 


.0 


0.99 + 0.00 


9.2 + 0.9 


14.1 + 1.2 


86 


27.8 + 0.5 


2.39 + 0.40 


Abelll942 


1 


.05 + 0.06 


4.3+ 1, 


.0 


0.44 + 0.01 


2.7 + 0.6 


62.2 + 25.9 


4 


9.3 + 1.4 


1.57 + 1.21 


Abell2104 


1 


.22 + 0.08 


6.1 + 1, 


.6 


0.68 + 0.14 


5.8 + 0.8 


38.1 + 9.0 


8 






Aboll2111 


1 


.07 + 0.10 


4.5+ 1, 


.5 


0.74 + 0.07 


7.3 + 2.5 


35.9 + 10.9 


129 


33.0 + 2.8 


3.22 + 2.43 


Abell2163 


1 


.38 + 0.11 


9.5 + 2, 


.5 


2.33 + 0.03 


12.0 + 1.2 


27.1 + 1.9 


160 


35.3 + 0.4 


3.76 + 0.37 


Abell2204 


1 


.34 + 0.07 


8.1 + 1, 


.6 


1.16 + 0.01 


8.7 + 0.6 


1.9 + 0.2 


< 3 


4.8 + 0.3 


< 0.1 


Abell2218 


1 


.14 + 0.08 


5.1 + 1, 


.4 


0.72 + 0.01 


4.3 + 0.6 


52.2 + 8.7 


60 


18.9 + 1.0 


1.28 + 0.53 


Abell2219 


1 


.35 + 0.07 


9.1 + 1, 


.9 


1.65 + 0.03 


7.1 + 0.9 


31.4 + 5.5 


8 






Abell2259 


1 


.05 + 0.09 


4.0+ 1, 


.2 


0.50 + 0.04 


4.1 + 0.9 


28.0 + 7.1 


78 


24.1 + 1.7 


1.18 + 0.95 


Abell2261 


1 


.52 + 0.05 


12.9+ 1, 


.6 


1.46 + 0.13 


6.6 + 1.0 


12.9 + 2.3 


< 2 


14.3 + 1.0 


0.39 + 0.21 


Abell2390 


1 


.33 + 0.06 


8.6+ 1, 


.5 


1.48 + 0.01 


11.0 + 0.9 


2.9 + 0.1 


4 


11.1 + 0.9 


1.24 + 0.17 


Abell2537 


1 


.22 + 0.05 


7.2+1, 


.1 


0.86 + 0.06 


5.9 + 0.9 


12.3 + 5.2 


17 


8.4 + 1.3 


0.99 + 0.74 


CL0024. 0+1652 


1 


.30 + 0.10 


9.8 + 2, 


.7 


0.45 + 0.08 


3.1 + 4.7 


16.6 + 7.7 


24 


73.5 + 11.5 


6.46 + 5.23 


MACS J0717. 5+3745 


1 


.46 + 0.07 


16.6 + 3, 


.4 


2.35 + 0.03 


12.3 + 1.9 


36.6 + 7.9 


224 


23.9 + 0.9 


23.09 + 3.24 


MACS J0913. 7+4056 


0, 


.95 + 0.07 


4.0+ 1, 


.3 


0.53 + 0.02 


4.8 + 0.7 


2.7 + 0.3 


4 


4.3 + 1.0 


0.40 + 0.19 


MS0015. 9+1609 


1 


.60 + 0.06 


21.9 + 3, 


.2 


2.01 + 0.06 


13.4 + 1.9 


16.7 + 2.9 


41 


8.6 + 1.1 


0.58 + 0.38 


MS0440.5+0204 


0, 


.85 + 0.06 


2.2 + 0, 


.7 


0.24 + 0.05 


2.8 + 0.5 


11.5 + 3.4 


< 3 


19.6 + 4.2 


1.38 + 1.16 


MS0451.6-0305 


0, 


.95 + 0.10 


4.5+ 1, 


.7 


1.03 + 0.02 


7.8 + 1.0 


28.0 + 5.8 


28 


11.9+ 1.1 


1.44 + 0.82 


MS0906.5+1110 


1, 


.36 + 0.09 


8.7+ 1, 


.9 


0.87 + 0.03 


3.5 + 0.5 


27.8 + 8.1 


3 


17.1 + 1.2 


0.20 + 0.14 


MS1008.1-1224 


1, 


.06 + 0.05 


4.8 + 0, 


.9 


0.58 + 0.04 


7.3 + 3.1 


16.8 + 7.9 


10 


55.8 + 2.2 


4.17 + 2.33 


MS1231.3+1542 


0, 


.54 + 0.11 


0.6 + 0, 


.4 


0.14 + 0.00 


1.4 + 0.1 


31.4 + 4.4 


72 


6.9+ 1.4 


5.08 + 3.82 


MS1358.1+6245 


1, 


.12 + 0.09 


5.9+ 1, 


.6 


0.67 + 0.07 


7.6 + 0.9 


6.4 + 2.0 


4 


8.6+ 1.2 


0.34 + 0.29 


MS1455.0+2232 


1, 


.04 + 0.05 


4.2 + 0, 


.8 


0.56 + 0.01 


3.1 + 0.2 


4.7 + 0.2 


3 


4.9 + 0.2 


0.13 + 0.06 


MS1512.4+3647 


0, 


.85 + 0.18 


2.6+ 1, 


.8 


0.34 + 0.03 


2.1 + 0.7 


8.3 + 3.2 


6 


6.7+ 1.3 


1.30 + 1.09 


MS1621. 5+2640 


1, 


.19 + 0.07 


7.7+ 1, 


.8 


0.83 + 0.03 


5.4 + 0.8 


24.5+ 15.7 


41 


19.0 + 4.3 


7.47 + 5.43 


RXJ1347.5-1145 


1, 


.25 + 0.12 


9.3 + 2, 


.9 


1.63 + 0.01 


13.1 + 1.8 


2.4 + 0.2 


< 4 


12.6+ 1.4 


1.30 + 0.41 


RXJ1524.6+0957 


0, 


.87 + 0.12 


3.4+ 1, 


.8 


0.41 + 0.04 


2.7 + 0.4 


27.2 + 28.9 


22 


63.2 + 5.6 


22.92+ 15.12 



Note. — All quantities are measured at r^Q, except for P3/P0 power ratio, which is measured at r^ik, and CbcGi which 
is in Mpc. Mx is the X-ray hydrostati c mass. Kg is the cen tral entropy, j^bcg the X-ray peak to BCG offset, wx is the 
centroid shift. Lensing masses are from |Hoekstra et al.| ||2012|l . 



CALDB for the remainder of this paper, stressing that 
any changes to our results would be well within the sta- 
tistical errors presented were we to switch to a different 
calibration release. 

To be able to combine Chandra and XMM-Newton 
data, one must first ensure that they are consistent. We 
find that the following simple cross-calibration prescrip- 
tion is able to bring the data into self-consistency: 

^co™ctcd(^) - AcxoiE) (^^^ ^ (7) 

where ^ = gives the unmodified CALDB area, and 
C > has the effect of down-weighting the high energy 



effective area of Chandra. We find that setting ( = 0.07 
brings Chandra and XMM-Newton masses into agree- 
ment as shown in Figure [2] In either case, the intrinsic 
scatter between Chandra and XMM-Newton mass mea- 
surements at these fixed radii is certainly less than 10%, 
though inconsistent with zero at the 2a level. 

We chose to modify the Chandra effective area, and not 
the XMM-Newton effective area, based on the fact that 
the XMM-Newton has exhibited the least variation over 
the years, whereas Chandra has enacted larger 10-15% 
changes in its effective area calibration historically. We 
note that had we modified the XMM-Newton effective 
area to match that of Chandra, then we would have found 
in what follows that clusters no longer exhibit self-similar 
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Fig. 3. — Bimodality in the joint distribution of BCG offset and central entropy; contours show Hnes of constant probabihty density after 
the points have been smoothed with a 0.25 dex Gaussian. The top and right axes show the ID probability density for central entropy and 
BCG offset. Blue triangles show cool-core clusters and red triangles show non-cool-core clusters. The horizontal thin line shows our chosen 
division between cool-core and non-cool-core systems, while the vertical line shows our chosen division between low BCG offset and high 
BCG offset systems. 



behavior and that (a) those with obvious substructure 
would be the ones whose masses calculated assuming hy- 
drostatic equilibrium would agree with their weak lens- 
ing masses, and (b) that clusters with cool cores would 
have hydrostatic masses greater than their weak lensing 
masses. This uncertainty in the telescopes' effective areas 
must be viewed as a fundamental systematic limitation 
of X-ray astronomy at least as related to cluster science. 

2.8. Online Data and Regression Tool 

All data and analysis software used for this paper are 
available online at http : / / sf star . sf su. edu/cccp. Fits 
of scaling relations (i.e., the modeling of linear or power 
law relationships among measured quantities) are com- 
plicated by the fact that error in both coordinates makes 
ordinary analysis invalid. A detailed treatise of recent 
developments in the theory behind modeling 2D data 
with e rrors in both coordinates appears in Hog g et al.| 
(2010). These techniques allow the simultaneous estima- 



tion of slope, intercept, and intrinsic s catter in such rela- 
tions. We implement the methods of Hogg et al. (20101 
at the data website for this article. 

3. MEASURES OF NON-RELAXED STATUS 

The gas in all clusters of galaxies exhibits some degree 
of deviation from an idealized smooth, triaxial distribu- 
tion. Such deviation could come in terms of subclump- 
ing, asymmetry, or both. Its presence gives some clue 
as to the nature of its evolutionary history; for example, 



asymmetry could indicate either the beginning or the end 
of a merger event; subclumps could either be recently ac- 
creted small groups of galaxies, or surviving cold cores 
from recent mergers. 

Despite this ambiguity, objective measures of substruc- 
ture are helpful in arriving at quantitative estimates 
of departure from equilibrium. To begin, we employ 
two common and well-tested measures of substructure: 
power ratios and centroid shift variance. Power ratios 
are Fourier-space estimators of fluctuations in the overall 
cluster surface brightness distribution, while the centroid 
shift is a measure of the variance of the distance between 
the X-ray surface brightness peak (which is always well 
defined) and the centroid (which in a non-relaxed clus- 
ter often varies significantly as a function isophote used 
for its esti mation) . We refer the reader to|Buote & Tsai 



(1995); Po ole et al., (2006); Jeltema et al 



I ( |!j005[ [2ml 

he calculation 



and |Bohringer et al.| ( |2010p for details on ti 
of these estimators. 

As further tracers of the relaxed or nonrelaxed state of 
a system, we also consider the somewhat more straight- 
forward measures, central entropy and the X-ray to op- 
tical center offset. Low central entropies indicate the 
presence of a cool core, which tend to be associated (non- 
exclusively) with relaxed clusters. We define the central 
entropy as: 



Kq = X(20kpc) 



(8) 



In other words, the central entropy is defined as the de- 
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Fig. 4. — Correlation of four different substructure measures (central entropy Kg, BCG offset DbcG: X-ray centroid variance wx, and 
P3/P0 power ratio) against each other. Blue triangles show cool-core clusters and red triangles show non-cool-corc clusters; blue circles 
show low-BCG-offset systems and red circles show high-BCG-offset systems. 



projected entropy profile evaluated at a radius of 20kpc 
from the cluster center. 

Similarly, the distance between the brightest cluster 
galaxy (BCG) and the X-ray surface brightness peak can 
be a good predictor of relaxed state, with lar ge shifts in- 



dicat ing ongoing or residual merger activity ( Poole et al. 
20071. We measure this distance via simple astrometry 



on X-ray and optical images, and call it Dbcg- 

One would expect relaxed halos to be more represen- 
tative of idealized halo growth models. Hence we expect 
scaling relations among the various thermodynamic and 
dark matter parameters to be tighter for clusters selected 
on the basis of the more well-behaved substructure indi- 
cators. We also expect the most powerful substructure 
measures to be correlated with each other. 



3.1. Correlations among measures of substructure 

We explore the possibility of whether our substructure 
measures show inherent correlation. The presence of such 
correlations, particularly when involving both X-ray and 
optical data, can serve as road maps towards our goal of 
quantifying departures from equilibrium as economically 
as possible. We use the Spearman's rank correlation co- 
efficient, with bootstrap resampling for determining Icr 
uncertainties. 

The relationship between central entropy and BCG off- 
set is the most significant correlation in our sample. This 
also happens to be the most interesting correlation due 
to the relative ease of deriving central entropy and BCG 
offset from observables. Figure |3] shows that the two sub- 
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structures measures appear to form a two-peaked joint 
distribution, with low central entropy, low BCG offsets 
in one corner, and high central entropy, high BCG offset 
clusters in another. The dividing line is best seen as a 
curve with equation 



Ko = 7keVcm^ 



Mpc 



-1/2 



(9) 



The high correlation coefficient between Kq and -Dbcg 
appears to be due to bimodality: when we calculate the 
correlation coefficient separately for either cloud, we find 
that the clouds individually do not contain significant 
internal correlation. Though the above formula offers 
the most clean separation between the two clouds, most 
of the separation can be captured by imposing cuts in 
entropy, or, somewhat less cleanly, in BCG offset. 

For this reason, throughout the rest of the paper, 
we introduce a labeling system that represents cuts in 
these two most easily measured substructure estimators. 
We use blue triangles to indicate Kq < 70 keV cm^ 
("cool core systems" or CC), and red triangles to in- 
dicate Kq > 70 keV cm^ ("non-cool-core systems" or 
NCC). This nomenclature is based on the fact that of 
70 keV cm^ corresponds to a cooling time of « 1.5 Gyr; 
most cool core clusters have central cooling times below 
this value. 

Similarly, we use blue circles to indicate systems with 
Dbcg < 0.01 Mpc ("low BCG offset systems") and 
Dbcg > 0.01 Mpc ("high BCG offset systems"). 

In Figure |4j we look for inherent correlations among 
the other various indicators of substructure. Strong cor- 
relations exist between the BCG offset Dbcg , the central 
entropy Kq, the X-ray centroid shift wx at t^q, and 



the P3/P0 ratio at 



2500 



(in measuring the latter two, 



we cut out the central 0.15 to avoid dilution of the 
signal by the cool core). Interestingly, the P3/P0 ratio 
measured at r'^Q (instead of t^qq) showed much larger 
scatter (presumably due to noise) and proved much less 
tightly correlated with the other substructure measures 
than the P3/P0 ratio at r^Jf,. 

In particular, it should be noted that P3/P0 exhibits 
almost as strong a correlation with BCG offset as does 
central entropy, though there is no evidence for bimodal- 
ity. For non-cool-core clusters, the P3/P0 is significantly 
more correlated with BCG offset than is the central en- 
tropy. This is quite a surprising result, since P3/P0 
traces cluster dynamics outside the cool core, whereas 
the central entropy is more sensitive to the inner parts. 

The BCG correlation trends are consistent with the 
well-known tendency of cool cores to occur in smoother 
(i.e. more relaxed, hence lower wx, low power ratio) clus- 
ters where a BCG sits clos e to the bottom of the potential 
well ( [Bildfell et"al]|2008[ ). This demonstrates the tight 
quantitative Imk between these completely independent 
X-ray and optical indicators of substructure. 

4. THE Lx-Tx RELATION 



Similarly to previous studies (e.g. Morandi et al.|2007 



Pratt et al. 2010 Mittal et al 



luminosity-temperature (Lx — J 



2011[ ), we find that the 
'x ) relationship exhibits 



a significant scatter of ~ 50% when the core of the cluster 
is included — a scatter which is diminished considerably, 
to 36%, when the core is excised. This effect is due to 



the overall non-self-similarity of cluster cool cor es in com 
parison to t he regions outside the cool core (e.g Vikhlinin 
et al.||20d6 ). When the core is not excised, the cool-core 
clusters lie significantly above the non -cool-core clusters, 
an effect first noted by Fabian' (1994 ) and subsequently 
studied in detail by |McCarthy et al.^ (2004, ). 

In Figure [5] and Table [5] we show tiiat when we include 
all cluster emission, the residuals of the Lx — Tx relation 
show a strong and significant correlation with both the 
central entropy of the cluster and the centroid shift wx 
(we choose because of the four measures discussed in 
^3.1 it offers the strongest correlation). However, when 



we cut out the central 0.15 r^Q , the distinction disap- 
pears, and the cool-core and non-cool-core clusters be- 
come indistinguishable in terms of entropy as well as wx ■ 
This is an example of "irreversible scatter" — in other 
words, outside their cores, the clusters of galaxies in our 
sample have "forgotten" the cause of the intrinsic scatter 
in the Lx-Tx relation, or at least sufficiently so that 
the substructure measures cannot be used to tighten the 
relations any further. 

5. LENSING MASS-OBSERVABLE RELATION 

The mass-observable relationship is an important in- 
gredient in the determination of the cosmological param- 
eters with clusters of galaxies. Because the mass func- 
tion is the ultimate connector between the cosmological 
parameters and the data, finding accurate mass proxies 
using multiple methods and wavelength regimes is im- 
portant. Comparison of X-ray derived observables with 
weak gravitational lensing masses, which do not require 
the assumption of hydrostatic equilibrium, has proved a 
fruitful path towards this end (e.g. ,Mahdavi et al.,|2008[ 



Okabe et al. 



2010 Jee et al. 2011 . We hst our results 



for several different mass-observable relations in Table [Sj 

5.1. Temperature, Gas Mass, and Pseuo-Pressure 

We begin by examining the lensing mass-gas tempera- 
ture relationship in Figure [6} while exhibit i ng significant 
intrinsic scatter (jVentimigUa et aL 2008 Zhang et al 
J iMantz et al ' 
n- 



2008 



WU), the M-T relation is stiU a 



worthwhile keystone for comparison with previous work. 
We find that the relationship is consistent with being self- 
similar, with a larger scatter and uncertainty at lensing 
'"soo than at X-ray r5oo. Regardless of whether we con- 
sider the cool-core or the non-cool-core subsamples, the 
scatter is roughly 46%. The scatter drops dramatically 
to 17% ± 8% when we use X-ray rsoo because of the in- 
herent correlation between the gas temperature and X- 
ray ^500 itself, which we do not attempt to model. The 
normalization derived for the mass-temperatur e relation 
is consistent wit h previous wo r k, for exa mple Pedersen 
& Dahlc (2007j), [Henry et al] ( |2009[ ) and [Okabe et al 
(2010|.^ 

Table [H also shows similar results for the core-excised 
X-ray luminosity-lensing mass {Lx — J^wl) relation. 
The instrinsic scatter (35%±13%) is consistent with that 
of the mass-temperature relation, and as before, the scat- 

than at , again likely 



ter is dramatically lower at 



500 



'500 



and Lx which 



due to internal correlation between ri, 
we do not model. 

Far more impressive is the gas mass-lensing mass re- 
lationship. The gas mass has been shown in previous 
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work to be a use ful mass proxy ( |Mantz et a"L]|2010| [OE 



abe et al 
clusters 



J] 

of 



2010) 



■essentially, the assumption that rich 
galaxies have the same gas fraction is turn- 
ing out to be a re markably robu s t one . We improve the 

significance of the Okabe et al. (2010) finding with our 

~WT — 



sample of 50 clusters: at 



' 500 



the gas mass is consistent 



with being proportional to the lensing mass, with a log 
slope of 1.04 ± 0.1, and a normalization implying a gas 



and its residuals compared to centroid shift variance wx ■ 
(bottom panels) same as top, except that the inner 0.15 has been removed. The residuals are uncorrelated with all four substructure 

measures. Blue triangles show cool-core clusters and red triangles show non-cool-core clusters. 

tity (lensing mass) is measured using optical data and 
the other quantity (gas mass) is measured using X-ray 
data. The issue is that the aperture itself, r^Q, depends 
on the lensing mass, and therefore, by choosing the same 
aperture for the gas mass, we might introduce a correla- 
tion that produces ar tificially low scatter. This effect was 
described in detail by Becker & Kravtsov (2011 ) who find 
that such correlations can result in the measured scatter 
being w 50% smaller than the true scatter. 

However, using a physical aperture of 1 Mpc com- 
pletely takes away any possibility of covariance between 
the two axes. In Figure [Sj we truly have two statistically 
independent observations, and yet the intrinsic scatter 
remains remarkably low, 16% ± 7%. The fact that the 
scatter does not change when switching to a fixed phys- 
ical aperture is reassuring. The la scatter uncertainties 
are just large enough to accommodate the scatte r under- 
estimate predicted by Becker & Kravtsov (2011) (e.g. if 
the "true" scatter at both and 1 Mpc is 20%, our Icr 
errors would be consistent with a 50% scatter underesti- 



fraction /, 



gas 



0.12 ±0.01. 



We find a low scatter of 15 ± 



for the M, 



gas 



Mr 



relation (Figure [7j) for all clusters, regardless of dynami- 
cal state. Interestmgly, the same scatter holds regardless 
of whether we use lensing or a fixed aperture of 1 
Mpc. 

This low scatter at fixed radius is important. Recently, 
sophisticated treatments of the covariance between the 
axes i n the mass-observ able relation have become pos- 
sible (Hogg et al. 20101. Specifically, in the case of gas 
mass and lensing mass measured at , there is a subtle 
correlation between the two axes, even though one quan- 
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Fig. 6. — The mass-temperature relationship at lensing rsoo (left) and at X-ray rsoo (right). The latter shows less scatter due to the 
intrinsic correlation of X-ray rsoo with temperature. Blue triangles show cool-core clusters and red triangles show non-cool-core clusters. 




Fig. 7. — The gas mass-lensing mass relationship at lensing rsoo- Blue circles show low-BCG-offset systems and red circles show high- 
BCG-offset systems. Most of the low BCG offset systems are also low central entropy clusters. 
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Fig. 8. — Lensing mass vs. gas mass at a fixed physical radius of 
1 Mpc. Blue circles show low-BCG-offset clusters and red circles 
show high-BCG-offset clusters. The relation retains the low scatter 
of the relations at fixed density contrast. 



;i6g slope: 6.56±'0,'08 
I Scatter: 23%±6% 




scatter we can expect due to geometric errors from the 
assumption of spherical geometry. Thus deviations from 
spherical symmetry can produce scatter we observe in 
the cool-core Afgas — M]^ relation, and as a result, we can 
begin to claim that we are approaching a full accounting 
of all sources of systematic error in the mass-observable 
scaling relation. 

We note that the BCG offset works as well as central 
entropy in identifying the low-scatter subsample. This 
is an interesting result, because of our four substructure 
measures, BCG offset is by far the least expensive to 
calculate, in that it does not require X-ray temperature 
(spectral) information — a set of X-ray and optical im- 
ages is sufficient to calculate Dbcg ■ However, it is worth 
noting that while the low BCG offset and cool-core sub- 
samples have significant overlap, they are not precisely 
the same, and the two cuts trace two different types of 
equilibrium (dynamical and thermal, respectively). 

Another frequently used mass proxy is Yx , the pseudo- 
integr ated pressure first pioneered by |Vikhlinin et al" 
( 2006 1 ; being the product of the gas mass and the core- 
cut temperature at r^Q, Yx is directly comparable to 
the integrated Sunyaey-Zel'dovich compton Y parame- 
ter ( [Plagge et aL]|2010[ [Andersson et alj[2011[ ) . 



We show the Yx-Ml relation at for our sample 
in Figure [9] we find consistency with the expected self- 
similar slope of 0.6, but slightly higher intrinsic scatter 
to the gas mass when used as a mass proxy: the overall 
intrinsic deviation is ~ 23% ± 6% regardless of whether 
we use the entire sample or the cool-core subsample. 

One might be tempted to argue that gas mass is a su- 
perior mass proxy to Yx, not simply because of its ease 
of calculation and comparable overall intrisic scatter, but 
also because of the systematically lower intrinsic scatter 
that comes about when only cool-core clusters are con- 
sidered. However, this discrimination between relaxed 
and non-relaxed clusters is perhaps not optimal in a cos- 
mological context, where uniformity of scatter across the 
entire sample is important. Where uniformity is most 
important, Yx is a superior choice to gas mass, because 
as we show in Table [3] it has uniform scatter regardless 
of cluster central entropy or BCG offset. 

Finally, it is instructive to compare Yx with its ra- 
dio counterpart, the cylindrical ly integrated Sunyaev- 



Hoekstra et al. (20121 



Y f<r'"- 



Zel'dovich (SZ) pressure Ysz 

consider direct correlations between Vsz from the Planck 
mission and projected weak lensing masses; they find an 
intrinsic scatter of 12 ±5% at projected 7'25oo- As a point 
of comparison, when we conduct a similar exercise on 
spherically determined Yx and Mwl (both measured at 
mate at and no scatter underestimate at 1 Mpc). In spherical r^po), w e find an intrinsic^scat ter^of 18%±6% 
Table [3] we also list the performance of Yx, Lx, and Tx, 
measured at fixed physical radius of 1 Mpc, as predic- 
tors of Mvkl(< IMpc). Overall, we find little difference 
between the intrinsic scatter at 1 Mpc compared to 



Fig. 9. — Lensing mass vs. pseudo-pressure Yx- Blue circles 
show low-BCG-offset clusters and red circles show high-BCG-offset 
clusters. 



consistent with the Hoekstra et al. ( 2012 ) SZ comparison 



500 



5.2. 



Regularity of cool core and low BCG offset clusters 

Another point of particular importance is the fact that 
for the cool-core clusters, the Icr scatter is < 10% (the 
scatter is < 6% if we cut on BCG offset instead) — 
these numbers are low enough to be cons istent with zero. 
Simulations and analytical work (e.g|Becker fc Kravtsov 



5.3. Predicting with fixed aperture mass proxies 

for surveys 

In a blind X-ray survey, the aperture 



500 



or even 

I may not be easily available. For cosmology, we still 
need to know M^qq. It is therefore useful to investigate 
whether one can directly predict M^^^ without the need 
to calculate overdensity radii ta for the various X-ray 
observables. For example, a wide-field all-sky X-ray 
survey may be able to measure hundreds of thousands 
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of gas masses within fixed physical apertures, but lack 
the photon statistics to allow for the calculation of X-ray 
overdensity ra dii. 

In Figures [TO] we consider this situation, plotting 
^wl(< '"50(f) against gas mass and Yx measured within 
a fixed radius of 1 Mpc. As expected, the slopes now 
deviate from self-similar, and the intrinsic scatter is con- 
siderably higher than in Figures [7] and [9] However, in- 
terestingly, Yx exhibits somewhat less scatter (29%) in 
this "mixed" scaling relation than does gas mass (40%). 
In surveys with poor photon statistics where no X-ray 
or weak lensing estimates of tsoq are readily available, 
Yx may constitute a better mass proxy. The results are 
summarized in Table [S] 

These data leave us with the perhaps dispiriting re- 
sult that low (< 10%) scatter X-ray mass proxies may 
be derived either at fixed physical radii, yielding total 
mass estimates within fixed physical radii; or they may 
be derived at fixed overdensity radii, yielding total mass 
estimates within fixed overdensities. But it seems diffi- 
cult to achieve very low scatter without either commit- 
ting to fixed physical radii (straightforward to measure, 
but more difficult to use for cosmology) ; or to fixed over- 
density radii (difficult to measure, but more useful for 
cosmology) in both axes. 

5.4. Lack of Correlation with Substructure Measures 

We have already argued that the intrinsic scatter in the 
lensing mass to X-ray observable relations is potentially 
fully accounted for by the triaxiality of the clusters; nev- 
ertheless, it is still useful to consider whether the scatter 
in such relation may be further minimized via correlation 
with measures of substructure, at least as an empirical 
means to gauge the effect of this triaxiality. However, 
we find that none of the substructure measures — BCG 
offset, central entropy, centroid shift variance, or power 
ratio — have any significant correlation with residuals in 
the mass-observable relation. We note that Marrone 



et al. (2012) did find a residual correlation with BCG 
elliptlcity m the relationship between weak lensing mass 
and the integrated Sunyaev- Zel'd ovich effect signal Ygz- 
We examine a similar relation in §6.2| 

While it is possible to minimize the scatter in the mass- 
observable relation by conducting a cut on central en- 
tropy, it is not possible to "correct" this scatter for the 
non-cool-core clusters by utilizing any of the four sub- 
structure quantifiers we consider or even BCG ellipticity. 

6. DEVIATIONS FROM HYDROSTATIC EQUILIBRIUM 
6.1. Hydrostatic Mass Underestimate 



Mahdavi et al. ( 2008 ) argued that a subsample of the 

that 



500 



clusters discussed here have X-ray masses at 
are on the average 15% lower than their lensing masses at 
r^Q. This discrepancy may be attributed to deviations 
from hydrostatic equilibrium due to residual gas motions 
and incomplete thermalization of the ICM; the fact that 
hydrostatic masses ten d to und erestimate the true mass 
was first discussed by Evrard ([ i990| ) and co i itinues to 



be importan t in simulations (Lau et al 



let aLH201 1| e.g. 



2009 ^ 

and observations of distan 



Cavaliere 
clusters 

Andcrsson et al. 2011 Jee et al.||2011|). While biases in 



biases are only ?^ 3 — 5% for our an alysis. (Becker & 
Kravtsov||2(nil jHoekstra et al.||2012| 



In l*'igure we extend our results to the full sample 
of 50 clusters. The larger size of the sample allows us to 
resolve differences between cool-core and non-cool core 
clusters. We find that cool-core clusters and non-cool- 
core clusters do not exhibit the same level of departure 
from hydrostatic equilibrium. 

Cool core clusters have hydrostatic masses that are 
proportional to their weak lensing masses at all radii. 
The AIx — M]^ relation for this subsample has a small 
scatter (< 20%), about the right level for all the scatter 
to be accounted for by triaxiality. Overall, we find that 
cool core clusters are consistent with having no difference 
between their X-ray and weak lensing masses. 

The picture is dramatically different for non-cool-core 
clusters. In these systems, we find a roughly constant 
hydrostatic mass to lensing mass ratio of 80%, regard- 



less of whether we look at 



'500 



'2500- 



Our results 



are consistent with N-body gas dynamical simulations 
as shown in Figure [12] and Table [4] Broadly, these re- 
sults are consistent tEe hydrostatic mass underestimates 
predicted by gasdynamical simulations that ac c ount for 
unthe r malize d gas , such as | Nagai et al. (2007), Jeltema 
et al. ( 2008 ) and |Lau et al 
non-cool-core clus 



We find that the 
;ers populate the lower end of the re- 



gion allowed by these simulations, whereas the cool-core 
clusters populate the region where X-ray and true mass 
agree within 10%. Of these simulations 'Jel tema et al. 
( [2008) ) is the most consistent with our measured 20% av- 
erage mass underestimate for disturbed systems. 

6.2. Correlation with BCG Ellipticity 

Finally, we consider the question of whether BCG el- 
lipticity is correlated with difference s in X-ray and w eak 
lensin g rn asses, as suggested by , Marrone et al. (2012). In 
figure 13 we show M^/M-wl at r25oo and r5oo, plotted 
against CFHT ellipticities measured at 30kpc. We find 
that cool-core systems are consistent with Mx/Ml = 1 
{x^/v — 18/14); whereas non-cool-core systems are 
definitively not consistent with M^/Ml = 1 {x^ /v = 
70/29). 

Therefore, for non-cool-core systems, we find a good 
correlation between BCG ellipticity and the X-ray to 
weak lensing ratio ratio at r^^jj^, and a weak correla- 
tio n at rsoo. Whi l e this is similar to the trend found 
by Marrone et al. (2012), there is a difference in that 
our cool-core systeins do not appear to participate in the 
correlation. Furtherm ore, also in apparent contrast with 
Marrone et al. (2012), our correlation becomes less sig- 



gravitational lensing masses could in principle also affect 
the X-ray to weak-lensing mass ratio, such systematic 



nificant at r^^ . We interpret this result as suggesting 
that while cluster orientation plays some role in low X- 
ray to weak lensing mass ratios, it is not the only agent 
at work in this complex relationship (indeed, the hydro- 
static mass underestimate must also play a role). 

We note that it is not altogether surprising that the 
trend of ellipticity with Mx/Ml for cool core clusters 
is insignificant. We have shown in ^ that our X-ray 
and weak-lensing mass es are consistent for this sub- 
population (in contrast, Marrone et al. (2012) contained 
several relaxed clusters with significant SZ to weak lens- 
ing mass discrepancies). 
Furthermore, it is difficult to untangle the effects of 
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Fig. 10. — The lensing mass at rsoo vs. gas mass and Yx measured at fixed radius of 1 Mpc. Blue circles show low BCG offset clusters, 
while red circles show high BCG offset clusters. 
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Fig. 11. — The relationship between hydrostatic mass and lensing mass at r|^^(^ ('e/^) and t^^^q (right). Blue triangles show cool-core 
clusters and red triangles show non-cool-core clusters. Cool core clusters tend to nave hydrostatic masses that agree with lensing masses; 
non-cool-core clusters tend to exhibit the hydrostatic mass underestimate. The solid line indicates the best fit; the long-dashed line indicates 
the line of equality; the short-dashed line corresponds to the cool-core clusters, and the dotted line corresponds to the non-cool-core clusters. 
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Fig. 12. — The X-ray to Weak-Lensing mass ratio as a function of density contrast for cool-core systems (blue triangles) and non-cool-core 
systems (red triangles). The error bars are not independent because the data within r25 00 als o contributes to the measurement at rsoo- 
The shaded region shows the range of X-ray cluster mass underestimate as determined by |Lau et al.| Ij2009| . 




Fig. 13. — The X-ray to Weak-Lensing mass ratio as a function of BCG ellipticity for 43 BCGs with measureable ellipticities at 30 
kpc. Shown are non-cool-core systems (red triangles) and cool-core systems (blue triangles). The correlation is significant only for the 
non-cool-core systems at ^2500; there is a marginal correlation at rsoo. Cool-core systems do not participate in the trend. 
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elongation along the line of sight (which would chiefly 
bias weak lensing masses high) and non-hydrostatic gas 
(which would chiefly bias the X-ray masses low). We 
also stress that the trend is altogether absent at rsoo- 
However, empirically, we can point out that the non-cool- 
core clusters with the highest ellipticities have consistent 
X- ray and weak lensing masses, something corroborated 
bylMarrone et all (|2012|). 



7. CONCLUSION 

We examine archival X-ray data on a sample of 50 clus- 
ters of galaxies; most of the clusters have Chandra data, 
while roughly half have XMM-Newton data of good qual- 
ity. All clusters have CFHT weak gravitational lensing 
data from either the CFH12k or the Megacam instru- 
ments. 

In attempting to combine Chandra and XMM-Newton 
data to maximize both effective area and spatial reso- 
lution, we confirm previously reported systematic cali- 
bration differences between the two observatories. Using 
multiple calibration releases, we find a 15% systematic 
difference in hydrostatic masses between Chandra and 
XMM-Newton. Reassuringly, there is no intrinsic scat- 
ter between the masses for the two observatories, indi- 
cating that the issue is merely a matter of overall gain 
calibration and not a more serious spatially dependent is- 
sue. We develop an effective area correction that revises 
Chandra masses downward into agreement with XMM- 
Newton masses. This correction is only valid for high 
temperature (> 5 keV) clusters such as ours; at lower 
temperatures, the two observatories are more consistent 
due to the abundant prominence of X-ray lines. 

Using the Lx — Tx relation, we find that our sample 
is consistent with being randomly drawn from the same 
parent population as samples with we ll understood selec- 
tion f unctions, suc h as HIFL UGS (Relprich & Bohringer 
^002] ) and MACS ( [Ebeling et al.,201(7f : 

We examine several measures of substructure, includ- 
ing central entropy, BCG to X-ray peak offset, centroid 
shift variance, and power ratios. There is a significant 
correlation among all the substructure measures. The 
most strikingly correlated quantities are the BCG to X- 
ray peak offset (in Mpc) and the central entropy mea- 
sured at a radius of 20 kpc. The hint of bimodality in 
the joint 2D distribution of the BCG offset and central 
entropy indicates a complex connection between the ther- 
mal and dynamical relaxation times of galaxy clusters. 

Gas mass is by far the most robust predictor of weak 
lensing mass, with < 10% scatter for cool-core clusters 
and 14%±6% scatter for the sample overall. It is followed 
by the X-ray pseudo-pressure, Yx, which has 23% ± 6 
intrinsic scatter for both cool core clusters and the sam- 
ple overall. The mass-temperature relationship has even 
higher scatter, 43%±21% for the sample overall. All scal- 
ing relations have slopes that are consistent with the ex- 
pected self-similar value. We also find that core-excised 
X-ray luminosity is somewhat better than temperature 
at predicting weak lensing mass, yielding 28% ± 18% in- 
trinsic scatter for relaxed systems. 

By comparing hydrostatic and weak gravitati onal lens- 
ing rnasses , we extend our earlier detection (Mahdavi 
et al. 2008) of non-hydrostatic gas, with associated de- 
viations from hydrostatic equilibrium, in X-ray clusters 
of galaxies. We are able to quantify the hydrostatic 



mass underestimate separately for cool-core and non- 
cool-core clusters. We find that cool-core clusters ex- 
hibit little or no difference between their weak lensing 
and X-ray masses; the hydrostatic mass underestimate 
is consistent with 0% at both r^^Q and at r'^Q. Non 
cool-core clusters, on the other hand, have fairly con- 
sistent, « 20% ± 10%, underestimates between the same 
radii. This is broadly consistent with N-body gasdynam- 
ical simulations of unthermalized gas. 

Except for the non-core-cut Lx-Tx relation, we do 
not find a significant correlation between the residuals 
in a given scaling relation and any of our four substruc- 
ture measures (central entropy, BCG offset, centroid shift 
variance or P3/P0 power ratio). We interpret this result 
as indicating that it is not possible to reduce the intrin- 
sic scatter in a scaling relation (other than the Lx — Tx 
relation) by applying corrections based on substructure 
measures to individual clusters. In essence, clusters of 
galaxies have "forgotten" the sources of their departures 
from self-similarity. This lack of correlation suggests that 
we may have accounted for most if not all the parameters 
that could affect the cluster selection function for cosmo- 
logical surveys, and that few if any "hidden" parameters 
remain. 

However, we do find a partial trend with cluster ellip- 
ticity: cool-core clusters have consistent X-ray and weak 
lensing masses at r^Q^] whereas non-cool-core clusters 

have increasing Mx{< ?-^TO)/^i(< ^2500) with BCG el- 
lipticity at 30 kpc. Clusters with low ellipticity BCGs 
are the most likely to have mismatched X-ray and weak 
lensing masses, while clusters with higher ellipticity are 
more likely to have concordant X-ray and weak lensing 
masses. We leave it to future studies to determine which 
combination of X-ray non-hydrostatic bias and lensing 
projection bias is contributing to this trend. 

We emphasize that the X-ray peak to BCG location 
offset is perhaps the most efficient among our inspected 
substructure measures. Selecting clusters based on low 
BCG offset is sufficient to guarantee scatter consistent 
with zero in the gas mass-lensing mass relation, at least 
for a sample as large or larger than ours. 

In summary, we find that cool-core clusters with Kq < 
70 keV cm^ or BCG offset < 0.01 Mpc are extremely 
well-behaved and regular systems with respect to their X- 
ray and lensing properties. However, it should be noted 
that the two cuts do not select the same subsamples, 
because low BCG offset is indicative of the dynamical 
equilibrium, whereas low central entropy is a result of 
thermal equilibrium. While there are clusters that are in 
both thermal and dynamical equilibrium, the overlap is 
not perfect. 

Clusters with > 70 keV cm^ show some intrigu- 
ing properties — such as tightly correlated P3/P0 power 
ratios and BCG offsets, a linear correlation between 
Mx/Ml and ellipticity, and consistently low X-ray to 
weak lensing mass ratios — but larger samples and more 
careful theoretical studies are required before we can 
learn how to use these relations to gain greater physi- 
cal insight into their evolution. 
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TABLE 3 

Mass Proxy Fits with Lognormal Intrinsic Scatter 





Proxy 


MwL 






Log Log 


Fra<;tional Scatter 


Proxy 


Aperture 


Aperture 


Sample 




Slope Intercept 


in Myvij at fixed proxy 






Relations at Fixed Overdensity 


in 


Proxy and Mass 




TJ"V8 keV 






all 


1 


97 ±0.89 1.04 ±0.06 


0.46 ±0.23 


TJ^VS keV 


^500 


^500 


all 


1 


42 ± 0.19 0.96 ± 0.02 


0.17 ±0.08 






rWL 
^500 


all 





54 ±0.12 0.81 ±0.04 


0.36 ±0.12 




'500 


'500 


all 





57 ±0.08 0.78 ±0.03 


0.27 ± 0.05 



A^Gas X '^(z) r^i^ r|^^ all 1.04 ±0.10 0.90 ± 0.02 0.15 ±0.06 

i^o < 70 kcV cm^ 0.91 ±0.20 0.89 ± 0.03 < 0.1 

Ko > 70 kcV cm^ 1.09 ±0.13 0.90 ± 0.02 0.18 ±0.09 

Dbcg < 0.01 Mpc 0.93 ± 0.13 0.89 ± 0.02 < 0.06 

^5bcg >0.01 Mpc 1.13 ± 0.18 0.90 ± 0.03 0.22 ±0.15 



YxEizf 



all 





56 


± 





07 





45 


± 





07 





23 


± 


0.06 


Ko < 70 kcV cm^ 





44 


± 





14 





53 


± 





12 





24 


± 


0.18 


Ko > 70 kcV cm^ 





62 


± 





09 





41 


± 





09 





21 


± 


0.08 


Dbcg < 0.01 





51 


± 





10 





49 


± 





08 





17 


± 


0.12 


Dbcg >0.01 





60 


± 





13 





42 


± 





11 





26 


± 


0.17 



Relations at Other Radii 



T^ys keV (keV) 

MGasE(z) 

T^"V8 kcV 
L'%^*E(z)-^ 

MoasEiz) 

Yx 



1 Mpc 



1 Mpc 


all 


1.10 ±0.57 


0.80 ±0.02 


0.15 ±0.11 




all 


0.26 ±0.07 


0.71 ±0.02 


0.19 ±0.04 


55 


all 


0.83 ± 0.14 


0.90 ± 0.03 


0.16 ±0.10 


55 


all 


0.40 ± 0.06 


0.48 ± 0.05 


0.12 ±0.04 


'500 


all 


3.04 ± 1.38 


1.03 ± 0.08 


0.46 ±0.31 


all 


0.58 ±0.15 


0.80 ± 0.04 


0.38 ±0.13 


all 


1.73 ±0.59 


1.20 ± 0.13 


0.39 ±0.18 


all 


0.80 ±0.15 


0.35 ±0.11 


0.28 ±0.14 



Note. — All proxies are fit against My/LEi^z) at the stated apertures. All masses are in units of IO^^Mq. The 
core-cut X-ray luminosity is in units of 10 erg s~ , and Yx is in units of 10" M0 keV. 



TABLE 4 

X-RAY TO Weak Lensing Mass Ratios 



Contrast Sample 


Mx/Ml 


Fractional Scatter 






in Mx at fixed M/, 



'500 



All 






0.92 ±0.05 


0.19 ±0.05 


Ko < 70 


kcV 


Clip 


1.11 ± 0.10 


< 10% 


Ko > 70 


kcV 


cni"^ 


0.19 ±0.06 




^^BCG < 


0.01 


Mpc 


1.04 ±0.07 


< 0.15 


Dbcg > 


0.01 


Mpc 


0.81 ± 0.07 


0.24 ± 0.07 


All 






0.89 ± 0.05 


0.20 ±0.05 


Ko < 70 


kcV 


9 

cm 


1.08 ±0.09 


< 9% 


Ko > 70 


keV 


cm"^ 


0.83 ± 0.06 


0.20 ± 0.06 


^bcg < 


0.01 


Mpc 


0.97 ± 0.07 


0.13 ±0.10 


Dbcg > 


0.01 


Mpc 


0.84 ± 0.06 


0.22 ± 0.07 


All 






0.88 ± 0.05 


0.21 ± 0.06 


Ko < 70 


keV 


cm^ 


0.97 ±0.10 


0.17 ±0.13 


Ko > 70 


keV 


cm^ 


0.83 ± 0.07 


0.22 ± 0.07 


^BGG < 


0.01 


Mpc 


0.85 ± 0.09 


0.22 ±0.11 


^BGG > 


0.01 


Mpc 


0.89 ± 0.07 


0.20 ± 0.08 
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